Coherent Destruction of Tunneling and Dark Floquet State 
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We study a system of three coherently coupled states, where one state is shifted periodically 
against the other two. We discover such a system possesses a dark Floquet state at zero quasi- 
energy and always with negligible population at the intermediate state. This dark Floquet state 
manifests itself dynamically in terms of the suppression of inter-state tunneling, a phenomenon 
known as coherent destruction of tunneling. We suggest to call it dark coherent destruction of 
tunneling (DCDT). At high frequency limit for the periodic driving, this Floquet state reduces to 
the well-known dark state widely used for STIRAP. Our results can be generalized to systems with 
more states and can be verified with easily implemented experiments within current technologies. 
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Two-state and three-state models are the simplest 
quantum systems. Despite their simplicity, they often 
provide very good approximations to describe realistic 
physical systems and are capable of revealing a variety 
of fascinating quantum effects. The understanding of 
their ubiquitous features are nowadays being exploited 
for the manipulation and control of quantum states of 
small systems involving single atoms, photons, or nano- 
devices Coherent destruction of tunneling (two- 

state models) and dark state (three-state models) are 
two of the elegant prototype examples where deep under- 
standing gained from quantum coherent effects in these 
simple systems are impacting the development of quan- 
tum technology in communication and computation. 

Coherent destruction of tunneling (CDT) was discov- 
ered in a periodically-driving double well system It 
describes a fascinating phenomenon whereby coherent 
tunneling between two wells (or the Rabi oscillation be- 
tween two states) is turned off by an externally enforced 
periodic level shift. Rs understanding is related to dy- 
namical localization [BJ, which occurs at isolated degen- 
erate points of the quasienergies [J, CDT has thus 
far generated significant interest, and has been theoret- 
ically extended into various forms ^ has also 
been observed experimentally in man y p hysical systems : 
including modulated optical coupler [19|], driven double- 
well potentials for single- particle tunneling [2fj| . and a 
single electron spin in diamond [2l|. More recently, it 
has also found application in tuning the tunneling pa- 
rameter of a Bosc-Einstein condensate [13, [23| . 

Dark state is often discussed in terms of a three-state 
system where two of them are coupled coherently to the 
intermediate state, as in the model system we study here. 
When all coupling fields are on resonance with their re- 
spective coupled pair of states, we can adopt the rotat- 
ing wave approximation and change into a suitable in- 



teraction picture with all coupling coefficients becoming 
time independent. In this case, there always exists a 
dark state, whose eigenenergy becomes uniformly zero, 
and the corresponding eigenvector contains no projec- 
tion onto the intermediate state. It is called dark as the 
intermediate state is an excited state capable of emitting 
photons. This type of dark state is also known as coher- 
ent population trapping [24[ , widely used in efficient pop- 
ulation transfer through the stimulated Raman adiabatic 
passage (STIRAP) protocol. It has become the theoret- 
ical basis for several well-established implementations of 
quantum control and rudimentary quantum information 
processing gates. 

In this Letter we report our surprising finding of an in- 
timate relationship between dark state and coherent de- 
struction of tunneling by studying a three-state system. 
In this system, two states are coherently coupled to an 
intermediate state and one of the two states shifts period- 
ically against the other two. We find that CDT also exists 
in this three-state system, where the dynamical tunnel- 
ing from one state to the other two is suppressed by the 
periodic driving over a range of parameters. However, 
this CDT for the three-state system has its own distinct 
feature: it is related to a dark Floquet state, which has 
zero quasi-energy and negligible population at the the 
intermediate state. Quite interestingly, this dark Flo- 
quet state reduces to the well-known dark state in a non- 
driving three-state A-systempij at high-frequency limit. 
Therefore, we call this CDT dark coherent destruction of 
tunneling (DCDT). These results can be generalized to 
A?"-state system. We also discuss a feasible experimen- 
tal scheme where the visualization of the DCDT can be 
achieved readily. 

Three-state system. The driving three-state system is 
described by the Schrodingcr equation (K = 1) 
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i— = -— sin(wt)c 3 + vc 2 , 

where c\, C2, and C3 are the amplitudes at three states 
|1), 1 2), and |3), respectively, v is the coupling con- 
stant between the neighboring states. Energy state |1) 
is shifted periodically against the other two with driving 
strength A and frequency oj. The normalization condi- 
tion E| =1 |cj| 2 = 1 is assumed. 

To investigate the tunneling dynamics, we solve the 
time-dependent Schrodingcr equation ([T]) numerically 
with initial state (1,0, 0) T . The evolution of the prob- 
ability distribution Pi = \c\\ 2 is presented in Fig. [T] for 
three typical driving conditions. For A/lu = (Fig. Q] 
(a)), we see that Pi oscillates between zero and one, 
demonstrating no suppression of tunneling. For A/uj = 
2.0 (Fig. [T] (b)), the oscillations of Pi are seen limited 
between 0.8 and 1, showing suppression of tunneling. At 
A/bj = 2.4 (Fig. [TJ (c)), Pi remains near unity, signal- 
ing a complete suppression of tunneling between energy 
states. This is the quantum phenomenon well known as 
CDT [|. 

We emphasize that what we find here is not a sim- 
ple re-discovery of CDT in a three-state system. The 
CDT in this three-state model has its own distinct fea- 
ture: the results shown in Fig. [TJ (b,c) indicate that the 
suppression of tunneling occurs over a wide range of sys- 
tem parameters. This is in stark contrast to the CDT in a 
two-state system [J] , which occurs only at isolated points 
of parameters. The widening of the suppression regime 
found in the driving three-state system is more clearly 
demonstrated in Fig. [1] (d), where the minimum value 
of Pi is used to measure the suppression of tunneling. 
When min(Pi) is not zero, the tunneling is suppressed as 
the population is not allowed to be fully transferred from 
state |1) to the other two states. It is clear from Fig. [1] 
(d) that the suppression occurs as long as A/uj ^ 0. For 
comparison, the results for the standard driven two-state 
system is plotted as dash dotted line in Fig. [T](d), where 
the extremely narrow peak width indicates that CDT oc- 
curs only at isolated points of parameters. 

There exists a fundamental reason why the CDT oc- 
curs at isolated system parameters in a two-state model, 
where the CDT is related to the degeneracy of quasi- 
energy levels [|| and the degeneracy usually happens 
only at isolated points. Therefore, the significant widen- 
ing of the suppression parameter range that we see in 
Fig. [TJ indicates that the CDT found here should have 
a different origin. To investigate this origin, we turn 
to the Floquet theory for a periodically-driving system. 
Similar to Bloch states for systems with spatially peri- 
odic potentials, the modulated system ([XJ) has Floquet 
states, (ci,C2,cs) T = (ci, £2, c^) T exp(— let), where e is 
the quasi-energy and the amplitudes (ci, C2, c~3) T are pe- 
riodic with modulation period T = 2tt/uj. 

Our numerical results of the quasi-energies and Flo- 
quet states for the modulated system fTJ) are plotted 
in Fig. [2J There are three Floquet states with quasi- 
encrgics ex, £2, and £3. We immediately notice from 




50 100 2 4 



t Mo 

FIG. 1: (color online) The evolution of the probability at 
state |1) Pi = |ci| 2 for the system (fTJ) for various driving 
conditions: (a) A/co = 0; (b) A/uj = 2.0; (c) A/u = 2.4. (d) 
The minimum value of Pi (solid line) as a function of driving 
parameter A/ui. The two-state results are plotted as a dash 
dotted line for comparison. The initial condition is {ci = 
1, C2 = 0, C3 = 0}. The other parameters are u = 10, v — 1. 

Fig. [UJa) that the quasi-energy £2 for the second Flo- 
quet state is always zero for all values of A/ to. We call 
this state dark Floquet state in analogy to the well-known 
dark state. This dark Floquet state stands out not only 
for its zero quasi-energy but also for its unique population 
distribution among energy states. We display the time- 
averaged population probability (Pj) == (j Q dt\cj\ 2 )/T 
for a given Floquet state (ci, C2, c^) T in Fig. [2Ib-d). The 
Floquet state with (|P/| 2 ) > 0.5 is generally regarded as 
a state localized at the j-th energy state. As seen in Fig. 
[2jc), the dark Floquet state has almost zero population 
at energy state |2) while the population at |1) (Pi) > 0.5 
. In other words, the dark Floquet state is localized at 
|1). The other two Floquet states have identical popula- 
tion distribution. Since all their populations (Pj) < 0.5, 
these two Floquet states are not localized. 

It is not difficult to see the suppression of tunneling 
seen in Fig.[TJis linked to the existence of the dark Floquet 
state. We expand the initial state in terms of the Floquet 
states 

(l,0,0) T = 6 1 |£ 1 )+6 2 |£2)+fe 3 |£3>- (2) 

During the dynamical evolution, the expansion coefficient 
hi evolve as bi exp (— ietf). In other words, |6,|s are time 
independent. We look at the case A/ui = 2.4, where 
|£2) has population one at state |1) while the other two 
states have zero population at |1). In this case, we have 
I&2I = 1 and bi = & 3 = 0, which corresponds to a com- 
plete suppression of tunneling from |1) to |2) and |3). 
For other values of A/ui, similar arguments can be made. 
This shows that the CDT observed in Fig. [1] has a differ- 
ent origin: it is the consequence of dark Floquet states. 
Therefore, we call it dark coherent destruction of tunnel- 
ing (DCDT). 
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FIG. 2: (color online) (a) Quasienergies versus A/u). Solid 
lines are for numerical results obtained from the original 
model (ft) and circles are for the approximation results given 
by the effective model © ■ Time-averaged populations for the 
Floquet state in the quasi-energy level (b) ei, (c) £2, and (d) 
£3. The other parameters are v = 1,uj = 10. 



Interestingly, this dark Floquet state can be reduced 
to the well known dark state in a non-driven three-state 
A-system at high frequency limit. Introducing the trans- 
formation c m = a m exp[±iAcos(ujt) / (2lo)] (+ for m = 1 
and — for m = 2,3) and averaging out high frequency 
terms, one can obtain a non-driven three-state system 

i ~^t = vJo{A/u)a 2 , 

i ~d~t = vJ o( A / UJ ) a ^+ va 3^ ( 3 ) 
.da 3 

l HT = m2 ' 

where Jq{A/w) is the zeroth order Bessel function. The 
famous dark state (also known as coherent trapped 
state) for Eq. ([3]) is given by (ai,02,a 3 ) T = 
-^(-v,0,vJ (A/uj)) t , where M = v 2 + [v J (A/uj)} 2 . 
This dark state corresponds to the dark Floquet state. 
Similarly, this dark state is always localized at state 
|1) as v > vJq(A/lo) and has zero population at state 
1 2). This state is completely localized at state |1) when 
Jo(A/u) = 0. We have computed the the eigenvalues of 
model ([3J and compared them (circles) with the quasi- 
energies (black solid lines) in Fig. H](a). The agreement 
is almost perfect. 

Generalization to N -state system. Our analysis above 
is given for a three-state system and the original CDT 
was found in a two-state system. These results can 
be generalized to an iV-state system, where one state 
is shifted periodically against all the other states. The 
equations of motion are 

dc 

= "(cj-l+Cj-H-O+^iWc,', (4) 



where Cj<a — c 3 >jv = 0. 
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FIG. 3: (color online) (a) The minimum value of Pi as a 
function of A/u) for N = 4 (dash dotted line) and N — 5 (solid 
line). The initial conditions are ci(0) = l,Cj(0) = 0(j 7^ 1). 
Quasi-energies versus A/u) for (b) N = 4 and (c) N = 5. Solid 
lines are for numerical results obtained from the model Q and 
circles are for the approximation results given by the effective 
model ([5]). (d) The time-averaged probability distribution of 
the Floquet state corresponding to e — in Fig. [3] (c). The 
other parameters are v = l,w = 10. 

The quantum dynamics of the driven TV-state sys- 
tems is investigated by direct integration of the time- 
dependent Schrodinger equation ((4]) with the state ini- 
tially prepared on the state |1). The CDT is found to 
exist. The minimum value of Pi = |ci| 2 as a function of 
Aju) is presented in Fig. [3] (a) for N = 4 and N = 5. 
When TV = 4, the CDT occurs at an isolated point of 
parameters (dash dotted line in Fig. [3] (a)), where two 
of the four quasi-energy levels become degenerate (Fig. 
[3] (b)). This is exactly the same as in the two-state sys- 
tem. When N = 5, the parameter range where CDT 
occurs is extended substantially (solid line in Fig. [3] (a)) 
as in the three-state model. Furthermore, this five-state 
system also has a dark Floquet state: as seen in Fig. [3] 
(c), one of the quasi-energies always equals to zero. This 
dark Floquet state has negligible population at all of even 
j-th states (Fig. [3J(d)). 

These numerical results with N = 4,5, together with 
the known results for TV — 2,3, clearly suggest that (i) 
the dark state and the associated DCDT exist in odd 
TV-state systems; (ii) the original CDT, which occurs at 
isolated parameter points, exists in all even-iV-state sys- 
tem. This general conclusion can be proved analytically 
at high frequency limit. 

Following the procedure used in the three- 
state system, we introduce the transforma- 
tion ci = a\ exp[— i J Asin(u)t)/2dt], Cj^\ = 
Oj-exp[i J Asm(u)t)/2dt], where aj(z) are slowly varying 
functions. Using the expansion exp[±iAcos(wt)/o;] = 
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^2k(±i) k Jk{A/uj) cxp(±iku>t) in terms of Bessel func- 
tions and neglecting all orders except k = for high 
frequency limit, we can reduce the coupled equations ((4]) 
to a non-driven model 

where a = (ai, ci2, cln) t ■ The matrix H is tridiagonal 

with H12 = H21 = V cS = vJ (A/uj), H n ,n+1 = Hn+l,n = 

v. The effective coupling constant v e s between state |1) 
and state |2) is tunable with the driving parameters. 

The eigenvalues and eigenvectors of the tridiagonal 
TV x TV matrix H enjoy some very interesting proper- 
ties, whose rigorous proofs can be found in Supplemental 
Material. 

(a) When TV is even, for any nonzero w e ff o- n d v, all 
the eigenvalues of the matrix H are nonzero while two of 
the eigenvalues are zero for v e s = 0. 

Remark: This means that when TV is even, two quasi- 
encrgy levels of the driven model Q are degenerate at 
isolated points where v g = vJo(A/uj) = 0. The CDT 
occurs at these isolated points. 

(b) When TV is odd, one and only one eigenvalue of H 
always equals to zero and, for the corresponding eigen- 
vector (wi, u>2, Wn) t of H, the inequality \w\\ 2 > 0.5 
holds for a finite range of parameters; for any other 
eigenvector {w\, W2, wjv) T of H , one has \wj\ 2 < 0.5. 

Remark: When TV is odd, the system always has one 
and only one dark Floquet state, which is localized over 
a finite range of parameters. Correspondingly, DCDT 
occurs over a finite range of parameters. 

Experimental observation. By mapping the temporal 
evolution of quantum systems into the spatial propaga- 
tions of light waves, the engineered waveguides have pro- 
vided an ideal platform to investigate a wide variety of 
coherent quantum effects [2^, [26| . The phenomenon of 
DCDT can also be observed with this kind of waveg- 
uide system. The discrete model (01 can be simulated 
by the light propagation in an array of TV waveguides 
placed closely and with equal spacing. Periodic driv- 
ing is realized by the harmonic modulation of the re- 
fractive index of the waveguides along the propagation 
direction [27], |28|. For our system, the periodic modula- 
tion of the first waveguide has a phase difference of it 
against the modulations for all other TV — 1 waveguides. 
When TV is odd, the DCDT can be readily observed with 
current experimental capacity [27l l28j . 

In summary, we find that the CDT also happens in a 
three-state quantum system, where one energy state is 
shifted periodically against the other two states. We call 
this type of CDT dark coherent destruction of tunneling 
(DCDT) as it is related to the existence of a dark Floquet 
state in the three-state system. The dark Floquet state 
has zero quasi-energy and negligible population at the in- 
termediate state. It reduces to the well known dark state 
of a non-driven three-state system. These results can be 
generalized to a periodically driven TV-state system. We 
have also pointed out that observation of DCDT is well 
within the capacity of current experiments. 
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Supplemental Material 

The TV x TV tridiagonal matrix H has the follow- 
ing non-zero matrix elements: H12 = H21 = v c g, 
H n , n +i = S n +i, n = v for n = 1, 2, ...,N - 1. v cS 
is tunable and v 7^ is fixed. The eigenvalues and 
eigenvectors of H have the following properties. 

Property 1. When TV is odd, one and only one eigen- 
value of H always equals to zero. 

Proof. Let Ai, A2, . . . , Xn be all the eigenvalues of matrix 
H, then Djq = det(-ff) = AiA 2 • • • Xn- It is easy to verify 
that Z?2 = — v 2 s , D\ = Z?3 = 0, and the relation Dn = 
—v 2 Dn-2{N > 3). Therefore, one has £>2fc-i = and 
D 2 k = (-l) k v 2k ~ 2 v 2 s (k = 1, 2, 3, . . .). When N is odd, 
Dn = A1A2 • • • Ayv = 0, which means that at least one 
eigenvalue equals to zero regardless the values of i/ e ff and 
v. For the zero eigenvalue, the eigen-equation is Hw = 
0, where w = (wi, W2, ■ ■ ■ , wn) t ■ When expanded, the 
equation is turned into the following equations: v e ffw 2 = 
0, v e gwi +VW3 = 0, vwj-i +vwj+i — (J = 3, 4, . . . , N — 
1), vwn-i — 0. There is only one non-trivial solution. For 
VcS 7^ 0, it is given by 

w=^=(w° 1 ,wl...,w° N ) T (6) 
VM 

with wi = (-ly-Wv/veB, w° 2k = 0, and w° 2k+1 = 
(-l)(N-2k-i)/2 whcrc k = i } 2,...,(N— l)/2 and M = 
X)j=i \ w j\ 2 i s t ne normalization factor. For v e g = 0, the 
solution isw=(l,0,...,0) T . This shows matrix H has 
one and only one eigenvalue equal to zero. 

□ 

Property 2. When N is even, for any nonzero v c f[, all 
the eigenvalues of the matrix H are nonzero while two of 
the eigenvalues are zero for v e s = 0. 

Proof. If i> ff 7^ 0, when TV is even, then Dn = 
(-l) N / 2 v N ~ 2 v 2 s = AiA 2 • • • A w ^ 0. Thus all the eigen- 
values of the matrix H are nonzero. When v c g = 0, it 
is obvious that Dn = A1A2 ■ ■ ■ An = 0. There are zero 
eigenvalues. With v s = 0, the tridiagonal matrix of H 



5 



is divided into two uncoupled subspaces, i.e., H = 0® F, 
where F is a tridiagonal matrix with nonzero elements 
F n , n +i = -Pn+l,n = v. It is clear that F possesses Prop- 
erty 1 and has only one zero eigenvalue. H thus has two 
zero eigenvalues. □ 

Property 3. If X is an eigenvalue of H with eigenvector 
(u>i, u> 2 , Wn) T , then—X is an eigenvalue of H with the 
corresponding eigenvector (w i 1 ,w i 2 t--,w i n ) t where w'j = 

(-iyw r 

Proof. The eigenvalue equation Hw = Aw can be written 
in the form Yi^ =l HijWj = Xwi, where Hij = when 
\i - j\ =0 and Ji - j\ > 2. Multiplying by (-l^A, 
we obtain Y, 1 j =1 Hij(—l) : >Wj = — X( — l) l Wi and have the 
proof. □ 

Property 4. When N is odd, for the eigenvector 



(wi,W2, ...,W]y) T of H belonging to X = 0, the inequality 
\u>i\ 2 > 1/2 holds for a finite range of parameters; For 
an eigenvector (wi, u>2, wn) t of H belonging to X ^ 0, 
one has that \u)j\ 2 < 1/2, whether N is odd or even. 

Proof. According to Eq. ((6]) , it is clear that the inequality 
|wi | 2 > 0.5 is valid only when (w/w ff) 2 > (N - l)/2. 
When A 7^ 0, the two eigenvalues A and —A are 
distinct and the corresponding eigenvectors are or- 
thogonal to each other. With Property 3, one has 
E& 1)/! Wi| a = E£r 15/ V*l a when N is odd; 

Efiil^-il 2 = X%H\w 2 k\ 2 when N is even. With the 
normalization condition E^L-Ju^l 2 = 1, we immediately 
obtain that \wj\ 2 < 0.5. □ 
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